Rotational properties of trapped bosons 
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Based on the Hellman-Feynman theorem it is shown that the average square radius of a cloud of 
interacting bosons in a parabolic well can be derived from their free energy. As an application, the 
temperature dependence of the moment of inertia of non-interacting bosons in a parabolic trap is 
, determined as a function of the number of bosons. Well below the critical condensation temperature, 

the Bose-Einstein statistics are found to substantially reduce the moment of inertia of this system, 
as compared to a gas of "distinguishable" particles in a parabolic well. 
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The study of a crowd of atoms in the same quantum state has gained great importance since the recent observation 
of the Bose-Einstein condensation of magnetically trapped gases of rubidium H, sodium and — albeit only with 
indirect evidence — lithium ||. According to Stringari 0], who draws a parallel with superfluid systems p], the 
rotational properties of a confined Bose gas provide a natural way to analyze the deviations from rigid motion due 
to condensation. He suggested that condensation effects of magnetically trapped gases should be observable in their 
moment of inertia, which he calculated approximately from the linear response properties of a Bose gas of harmonic 
oscillators without two-body interactions, relying on the semiclassical approximation for the excited states. 

Using the Hellman-Feynman theorem, we show in the present letter that the moment of inertia of a confined Bose 
gas can be derived directly from the dependence of its free energy on the strength of the confining parabolic potential, 
\Q ' even if two-body interactions are included. Indeed, for the following Hamiltonian 
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it is straightforward to show that the confinement contribution U con f = (Uconf) to the internal energy can be written 
as 

Uconf = ^Tr (e-^Hconf) = jN (cj 2 z (z 2 ) + lo 2 (x 2 + y 2 )) , (2) 



where Z is the partition function. Denoting the free energy by T , the Hellman-Feynman theorem states that 



relating thereby the average square radius (p 2 ) = (x 2 + y 2 ') of the system in the rry-plane (and hence its moment of 
inertia with respect to the z-axis) to partial derivatives of the free energy. This relation is valid both for distinguishable 
and for identical particles. If the free energy exhibits critical behavior, it will also be reflected in the average square 
radius of the Bose gas and consequently in its rotational properties. 

In order to estimate the moment of inertia quantitatively, we consider for simplicity an ideal cloud of non-interacting 
bosons trapped in a parabolic well characterized by a frequency w. With fi denoting the chemical potential and 
j3 = 1/kT where k is Boltzmann's constant, and T the temperature, the partition function in the grand-canonical 
ensemble becomes 
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where | (v + 1) (y + 2) is the number of different possibilities to occupy the v th level. This multiplicity depends on 
the symmetry and is considered here for the isotropic case to 2 = u> 2 . For cylindrical symmetry, e.g., it would be 1 
along the symmetry axis, and (y + 1) in the plane orthogonal to the symmetry axis. 

The expression (|4|) is the isotropic limit of the grand canonical partition function for a Bose gas of anisotropic 
oscillators ||, from which the subsequent derivation can be generalized to the anisotropic case to 2 ^ ui 2 . The 
subsequent discussion is limited to the isotropic Bose gas. The ground state occupancy and the thermodynamic 
properties including the aniso trop ic case have recently been analyzed in several papers J^— and are currently also 
the subject of some preprints 

The average number of bosons N — ^ ^- In S M and their internal energy U = — ^'^p + N/i in an isotropic parabolic 
confining potential are: 
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N = 22n u ; U = 22n u e„; n v = - (y + 1) (u + 2) - _ , (6) 

u=0 v=0 

where n v is the occupation number of the v til level. Furthermore the Hellman-Feynman theorem leads to the average 
square radius of the boson cloud 

-mNlr) = --_ — .lnS„ = — ^ > {v + l) (u + 2)e v - m r = (7 
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Consequently, the averaged potential energy becomes ^mw 2 N (^f 2 ^ = \\J, as could have been anticipated from the 
virial theorem. The moment of inertia relative to the "Z-axis is therefore given by: 

Q = mN ( X 2 +y 2 ) = \^ (8) 

which is the basic result of the present letter. Critical behavior of the internal energy (and consequently in the specific 
heat) of the model is directly reflected in the rotational properties of the boson gas. 

The evaluation of the internal energy for a given number N of particles, requires the knowledge of the chemical 
potential. By expressing the chemical potential as a function of the number of particles in the ground state 



the total number of particles N is obtained in terms of no: 
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However, as it turns out, this series is not very appropriate for numerical treatment. A much more efficient and 
numerically well-convergent expression results if each term in the denominator of (|l0|) is expanded in powers of 
e -f)hwv ^ w j 1 j c } 1 allows to perform the summation over v. The same procedure can be used for the internal energy. The 
results can be written in the following form, which is numerically stable: 

N = 52i=t {r^+Tj {i-e-^if = n o + J2e=i (^+i) {^Z^m^Tja - l) 5 

- 2 hW l^£=l ) (i_ e -3W)4 - 2 hW + 2 l^i=l {iZ+T ) (1 _ e -0fi^)<t J ■ 

In these expressions, I is the length of a cycle in the cyclic decomposition of the partition function obtained from 
the symmetrized density matrix for N oscillators without two-body interactions. 

The high temperature limit 8 — * of the internal energy and of the depletion of the condensate can now be analyzed 
directly. There is no need anymore to make the continuum approximation by introducing a parametrized form for 
the density of states ||. In this limit, and without making the continuum approximation, N — no becomes: 
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and similarly for the internal energy U. For no 3> 1, the summations yield 



fcr^i N-no = N(X-) + )+0{T) 
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This result for N ~ no has been derived (also without the continuum approximation) in and coincides with the 
result within an improved continuum approximation as derived in II and pj] |. The temperature expansion for the 
internal energy C7 is obtained in H using the continuum approximation. 

For arbitrary temperature, the ground state occupancy no and the internal energy U have to be evaluated numeri- 
cally from (|TT|) . This can efficiently be done by first solving for the quantity which is bracketed between and 1. 
The resulting ground state occupancy clearly exhibits the onset of a phase transition near the transition temperature 

Denoting by ©£> the moment of inertia of the "distinguishable" oscillators (i.e. without Bose-Einstein statistics), 
we also calculated the relative excess moment of inertia (O — Qd) /®d numerically. The results are shown in Fig. 1 
as a function of the reduced temperature for several values of N. For N of order 100 or more the effect of the Bose 
statistics on the moment of inertia is very pronounced, even reaching 75% for N — 1000 and 90% for N — 10000. 

The Hellman-Feynman theorem is also applicable to energy estimates based on the variational principle of quantum 
mechanics 16 17|], and it can be used for other models than the isotropic gas of bosons in a parabolic well as considered 




here. For instance, by applying the Hellman-Feynman theorem using the internal energy from the semiclassical 
approximation , one finds the same moment of inertia as was derived by Stringari Q] from the rotational response 
properties of this semiclassical model. Note that the functions Q (if) — J °° s v / (exp s — 1) ds, used in equation (9-10) 
and (13) of Ref. Q, are related to the Riemann £ functions: Q (2) = 2£ (3) and Q (3) = 6C (4) , and are a special case 
of the well-known Bose-Einstein functions g n (z) = Sfci z3 1 f 1 1 t ne properties of which are described in standard 
textbooks, e.g., in Ref. |ll|. 

It should be emphasized that the Hellman-Feynman theorem is valid in the canonical as well as in the grand 
canonical ensemble. Because the internal energy in the grand canonical ensemble with an average number of N 
bosons equals the internal energy in the canonical ensemble with a fixed number of N bosons fl9| , also the moments 
of inertia of both ensembles are the same. 

Having elucidated that the moment of inertia can be derived from the confinement energy in some models for 
the ideal or almost ideal gas, we did not discuss how quantitative results can be obtained in the presence of two- 
body interactions. Well aware of this limitation in the present letter, we like to point out that the inclusion of 
harmonic two-body interactions presents no conceptual difficulty. Their influence on the partition function and on the 
thermodynamical properties can be and has been obtained |l{| with a technique which we developed for treating the 
path integral for identical particles using the Feynman-Kac functional and invoking an appropriate linear combination 
of boson and fermion diffusion processes p0|-p2| . 

In summary, the partition function of an ideal Bose gas in a parabolic trap is has been obtained in the present 
letter, and the ground state occupancy no, the internal energy U and the moment of inertia have been derived. The 
rotational properties of the trapped gas are related to the internal energy through the Hellman-Feynman theorem. 
In the isotropic model we found the relation = |^ , implying that the critical behavior of the internal energy will 
manifest itself macroscopically in the moment of inertia. The semi-classical approximation [^0^1 turns out to provide 
the correct asymptotic limits of this exactly soluble model. Moreover, the path integral method allows to study 
this model as the zeroth order approximation in examining the role of the two-body interactions, using many-body 
diffusion and the Feynman-Kac functional. 
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Figure Captions 

Fig. 1 Relative excess moment of inertia (O — 0£>) /@d with respect to the z-axis as a function of T/T c for N = 10, 
100, 1000 and 10000. 
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